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Abstract 

We study a broad class of isotropic vacuum cosmologies in fourth- 
order gravity under the condition that the gravitational Lagrangian 
be scale-invariant or almost scale- invariant. The gravitational La- 
grangians considered will be of the form L = f(R) + k(G) where R 
and G are the Ricci and Gauss-Bonnet scalars respectively. Specifi- 
cally we take f(R) = aR 2n and k{G) = (3G n or k(G) = f3G In G with 
a, f3 and n arbitrary. We find solutions in closed form for a spatially 
flat Friedmann space-time and interpret their asymptotic early-time 
and late-time behaviour as well as their inflationary stages. 

One unique example which we discuss is the case of a very small 
negative value of the parameter f3 in the Lagrangian L = R 2 /3 + 
f3G In G which leads to the replacement of the exact de Sitter solution 
from L = R 2 (being a local attractor) to a power-law inflation exact 
solution a(t) = t p = t~ 3 ^ also representing a local attractor. This 
shows how one can modify the dynamics from de Sitter to power-law 
inflation by the addition of the G In G-term. 
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1 Introduction 



From the huge class of theories of gravitation which can be considered for 
describing and explaining the early evolution of the Universe, it is the sub- 
class of scale-invariant ones which plays as prominent role. The reason for 
this prominence is that almost all physical theories and their resulting cos- 
mologies have some limiting regime which is free from any scales. 

In the present paper we investigate a large class of vacuum cosmologies 
which include not only scale-invariant theories but also almost scale- invariant 
theories (we will define shortly the meaning of "almost scale-invariant theo- 
ries"). The Lagrangians which we consider are of the form 

L = f(R) + k(G), (1.1) 

where R is the curvature scalar and 

G = RijkiR 1 ^ 1 — ARijR 1 ^ + R 2 

is the Gauss-Bonnet scalar. By placing restrictions on the form of the func- 
tions f(R) and k(G) we will obtain theories which are scale- invariant and 
almost scale-invariant. We will focus on cosmo logical solutions to these scale- 
invariant and almost scale-invariant Lagrangians - in particular spatially flat 
Friedmann space-times. We obtain general features for these almost scale- 
invariant cosmologies, and for certain cases we are able to completely inte- 
grate the resulting Friedmann- like equations to confirm older results and to 
obtain new exact, closed form solutions in terms of the Friedmann metric 
scale factor, a(t). 

We first give the exact definition of what we mean by a scale-invariant 
gravitational Lagrangian. A theory of gravitation with a geometric 2 La- 
grangian L = L(gij) is called scale-invariant if a change of the scale in the 
metric (i.e. — > e 2c gij with c a constant) induces a change in the scale 
of the Lagrangian given by 

L(e 2c 9ij ) = e mc L{ 9l] ) , (1.2) 

where m is a constant. This has the consequence that if some particular 
metric gij is a vacuum solution then so too is the homothetically related 
metric e 2c gij. Einstein's theory of general relativity is scale- invariant in this 

2 We allow L to depend on the metric and its partial derivatives up to arbitrary order. 
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sense, but only if the cosmological term is absent 3 . Another example of a 
scale-invariant theory follows from the Lagrangian 

L = R 2n ■ f(G/R 2 ) 

with a constant n and an arbitrary (transcendental) function /. However, 
this class is still too large to be tractable. 

Next we define what we mean by almost scale-invariant. We call the 
Lagrangian L almost scale-invariant, if instead of eq. (1.2) we have the 
weaker condition 

L(e 2 %) = e mc L(g l3 ) + v(c) • H( 9ij ), (1.3) 

where m and c have the same meaning as above; v(c) is an arbitrary func- 
tion; H is a total divergence which therefore does not contribute to the field 
equations. Since H does not contribute to the equations of motion we find 
that the almost scale-invariant Lagrangians as defined in eq. (1.3) also have 
the feature that "if some particular metric is a vacuum solution then so 
too is the metric e 2c gij n . 

We now specialize eq. (1.1) to the case of Lagrangians which are scale- 
invariant. The following form of the Lagrangian which is a power law form 
in R and G with n, a and j3 constants 

L = aR 2n + (3G n . (1.4) 

is scale-invariant, and no further Lagrangians share this property. This La- 
grangian represents a closed set of functions i.e. all limits of elements of this 
class of functions are still within this class. 

We now expand our class of Lagrangians to include almost scale-invariant 
Lagrangians as defined in eq. (1.3). Since G is a total divergence, we can 
add some multiple of it to the scale-invariant Lagrangian eq. (1.4) like 

L = aR 2n + (3G n + 7 G = aR 2n + G ■ 0G n - 1 + 7). (1.5) 

The Lagrangian in eq. (1.5) is different from eq. (1.4) in that it does not 
represent a closed set of functions. This can be seen as follows: We take the 
limit n — > 1 in (1.5) while a remains constant and (3 = —7 = j3/(n— 1). We 
get 

lim aR 2n + (3G ■ — — — = aR 2 + (3G In G. 

n->l n — 1 

3 or is interpreted as part of the matter action 
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In arriving at this result we have defined n — 1 = e and G = e z and used 
lim — = lim — = lim = z = In G. 

n-*l n — 1 e^O £ e^O £ 

Result: The almost scale- invariant Lagrangians, as defined by eq. (1.3), 
expand the class of scale-invariant Lagrangian given in eq. (1.4) to include 
a Lagrangian of the following form 

L = aR 2 + (3GlnG + -fG. (1.6) 

Now we note that, since the argument of the logarithm in eq. (1.6) should 
be positive and dimensionless we should replace InG by ln(G/G ) where G 
is some constant G ^ 0. Further since G is a total divergence we can absorb 
the term 7G in eq. (1.6) into the GlnG-term via a suitable redefinition of 
Go- Thus the almost scale- invariant Lagrangian can be written as 

L = aR 2 + /3G\nG (1.7) 

with constants a and (3 fulfilling j3 7^ 0. The form of the Lagrangian given 
in eq. (1.7) that follows from the condition of being almost scale-invariant 
given in eq. (1.3) makes this case worth discussing since the condition of 
almost scale-invariant gives a very specific form of the Lagrangian. 

Before moving to the detailed calculations we give a very brief review of 
other work in this area. Cosmological models where the action depends on 
the Gauss-Bonnet scalar G are discussed in [1], [2], [3], and [4]. In [5], exact 
solutions for k(G) = G 13 are given. Further papers on this topic are [6], [7], 
[8], [9], [10], [11], [12], [13], [14], [15], [16], [17], and [18]. In [19], solutions for 
L = R+\f~G with a Friedmann scale factor of power law form (i.e. a(t) oc t p ) 
are given. In [20], the Lagrangian L = R n + (3G n ^ 2 is preferred. Further 
models are given in [21], whereas in [22] the case k(G) = G n + GlnG is 
discussed. 

Analogous models for /(_R)-gravity can be found in [23], where the case 
L = R 3 / 2 is related to Mach's principle. In [24] an exact solution for L = R 2 
is given. Further models are discussed in [25], [26], [27], [28], [29], [30], [31], 
[32], [33], [34], and [35]. In [36] the case f(R) = R 3 / 2 is studied. Newer 
models of this kind can be found in [37], [38], [39], [40], [41], [42], [43], [44], 
[45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [58], [59], 
[60], [61], [62], [63], [64], [65], [66], [67], [68], [69], [70], [71], [72], and [73]. 
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The conformal Weyl theory, especially the value of the perihelion advance 
in this theory, has been discussed in [74], [75], [76], [77], [78], [79], [80], [81], 
[82], and [83]. For theories in lower- dimensional space-times see e.g. [84], 
[85], and [86]. 

2 On the Gauss-Bonnet scalar 

The field equations for the Lagrangian L = f(R) are given by (see for example 
eq. (2.27) of [35] and the surrounding discussion) 

= L R R ij -g ij L/2 + g ij DL R -(L R )^ where L R = df/dR. (2.1) 

Since the case L = f(R) has been widely studied we will not go into further 
details here but simply refer the interested reader to the overview [35] . 

The case when L = k(G) is much less known than the case L = f(R) so 
we give some further details here. For the spatially flat Friedmann metric 
(given below in eq. (3.1)) the Gauss-Bonnet scalar G becomes 

G = 2Ah 2 (h 2 + h)= 24/i 4 (l + 7) , (2.2) 

where h is the Hubble parameter h = a/a and 7 = h/h 2 . For the Lagrangian 
k(G) with k G = dk/dG, the corresponding vacuum field equation is given in 
eq. (3.3) of [17] as 

= \g ij k{G) - 2k G RR ij + Ak G R l k R kj - 2k G R lklm R j klm 

-4k G R lkl3 R M + 2Rk% - 2g ij RUk G - ±R ik k%. k 
-4R jk k G . k + m ij Uk G + 4g ij R kl k G ., kl - 4R ikjl k G . kl . (2.3) 

Eqs. (A4), (A5) of [17], specialized to the space-time dimension n — 4, give 
rise to 

Rijki = Cijki + - {Rik9ji + RjWik — RiWjk — RjkQii) 

- g-R (9ik9ji - 9u9jk) ■ (2.4) 

Cijki is the Weyl tensor and we define 

U — U Uijki 
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and 

yij = R ikim R J kim + 2 R iklj R kl . (2.5) 

With this notation, the first (unnumbered) equation of the appendix of [17] 
reads 

C mm C 3klm = -5)C 2 . (2.6) 
Inserting eqs. (2.4) and (2.6) into eq. (2.5) we get 

yij = l g ij C 2 + l R 2 g ij _ RRij _ \ g i3 RklR ki + 2 R lk R\ . (2.7) 
4 6 2 

Further it holds that 

R m R ijkl = C 2 + 2R kl R kl - ^R 2 (2.8) 

and 

G = C 2 - 2R kl R kl + -R 2 . (2.9) 
With these notations we can rewrite eq. (2.3) as 

= -g^ (k(G) - Gk G ) + 2Rk% - 2g ij RUk G - 4R lk k' G . k 

- + 4R ij nk G + 4g^R kl k G;kl - AR ikjl k G . M . (2.10) 

Inserting k(G) = G n into eq. (2.10) we get with k G = nG'^ 1 

77 — 1 

= — g ij G n + 2nR(G n - l y tj - 2ng lJ RU(G n - 1 ) - AnR lk (G n ~ l )% 

-An (R jk (G n '% - R ij O(G n ~ 1 ) - g ij R kl (G 11 ' 1 ) . kl + R ikjl (G n ~ l ) . k ^j .(2.11) 

Inserting k(G) = G ■ lnG into eq. (2.10) we get with k G = 1 + lnG 

= --g ij G + 2R(\nG)' ij - 2g ij RU(\aG) - AR ik i\nG)% 
-AR jk (\nG)i + AR ij U(\nG) + 4g i] R k \\nG)- M - AR ikjl (\nG). kl . (2.12) 

3 General cosmological solutions 

In this section we use the background developed in the previous section to 
give a general study of spatially flat Friedmann space-times for the scale- 
invariant and almost scale-invariant theories cataloged in the introduction. 
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In subsection 3.1 we study the cosmology for the Lagrangian L = R 2n ; in 
subsection 3.2 we study the cosmology for the Lagrangians L = G n and L = 
G In G; in subsection 3.3 and 3.4 we study the cosmology for the combinations 
L = aR 2n + f3G n (scale-invariant) and L = aR 2 + (3G In G (almost scale- 
invariant) respectively. 

We start by setting up our system and notation. First, the cosmological 
metric we use is the spatially flat Friedmann space-time given as 

ds 2 = -dt 2 + a 2 {t) (dx 2 + dy 2 + dz 2 ) (3.1) 

with positive cosmic scale factor a(t). The dot denotes d/dt, h = a/a is the 
Hubble parameter, and R = Q(2h 2 + h) is the curvature scalar. 4 Without 
loss of generality we assume h > 0. If this is not the case then it is always 
possible to invert the time direction so as to get h > 0. If h appears in the 
denominator, this automatically includes the additional assumption, that 
h ^ 0. 5 It proves useful to define the function 

7 = h/h 2 = -± (I) (3.2) 



dt \h 

which shall be used to replace h in subsequent formulas. In terms of 7 we get 
R = Qh 2 (2 + 7). The deceleration parameter (i.e. q = — da/ (a) 2 ) is related 
to 7 via q = — 1 — 7. 

Sometimes it proves useful to use r = In a as an alternative time coordi- 
nate. With a dash denoting d/dr, we get with f = h the following formula: 

7'^ = ^-^ = ^ (3.3) 
1 dr dt dr h K J 

We now give some results which will be useful in dealing with the scale 
invariant Lagrangians of the form given in eq. (1.4) or almost scale- invariant 
Lagangians of the form given in eq. (1.7). First we note that, assuming 
a spatially flat Friedmann metric of the form given in eq. (3.1), that the 



4 If the metric signature is changed to (H ) then R must be replaced by —R, 

whereas G remains the same. In regard to the ambiguities in the sign of R often found in 
the literature we note that we define R such that for the standard 2-sphere we always get 
R> 0. 

5 This is not a real restriction, as a constant function a(t) is the trivial Minkowski space- 
time with h = 0, and solutions, where h{t) = at isolated points t only, can be matched 
by pieces with h 7^ 0. In other words: if h(t) = at isolated points t then these are always 
connected to regions where h{t) 7^ 0. 
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vacuum field equation for a Lagrangian of the form L = F(G,R), where F 
is a function of R and G, is (see eq. (15) of reference [5]) 

= GF G -F- 24h 3 F G + Q(h + h 2 )F R - QhF R (3.4) 

where F G = OF/dG and F R = OF /OR. Eq. (3.4) is the 00-component 
of the vacuum field equation for the Lagrangian L = F(G,R). All other 
components of the vacuum field equation are fulfilled if eq. (3.4) is valid. For 
the scale- invariant Lagrangian from eq. (1.4) 

L = aR 2n + (3G n 

eq. (3.4) becomes 

= a ■ T{2 + 7 ) 2n ~ 2 (2n( 7 2 + 3 7 + 2) - (2 + 7 ) 2 - 2n{2n - 1)( 7 ' + 4 7 + 2 7 2 )) 

+/3 • {n - l)2 n (l + 7 ) n ~ 2 ((1 + 7 ) 2 - n( 7 ' + 4 7 + 4 7 2 )) . (3.5) 

Before going in details we discuss some special cases. First we will exclude 
the case of Minkowski space-time so that we will have (with the possible 
exception of singular t-values) h > which implies that 7 is always well 
defined. Second, we exclude R — which gives 7 = —2 which in turn implies 
that a(t) oc t 1 ! 2 in eq. (3.1), i.e. one has a radiation-dominated Friedmann 
space-time, which is a self-similar geometry. Eq. (3.1) with a(t) oc t 1 ^ 2 is 
a solution for any L = R 2n with n > 1; is not a solution for < n < 1; 
is out of the range of applicability of the theory for n < 0. The remaining 
self-similar spatially flat Friedmann space-times for the metric in eq. (3.1) 
have a(t) oc t p with p ^ and p ^ 1/2. This scale factor a(t) = t v implies a 
constant 7 , namely 7 = —1/p. 

3.1 Cosmological solutions in i? 2n -gravity 

For L = R 2n , i.e., with a — 1 and f3 — 0, eq. (3.5) reduces to 

= -4 - 4 7 - 7 2 + 2n( 7 2 + 3 7 + 2) - 2n(2n - 1)( 7 ' + 4 7 + 2 7 2 ) . (3.6) 

For n = 0, no solution exists. The self-similar spatially flat Friedmann models 
are found by inserting 7 = — 1/p into eq. (3.6). This leads to 

= Ap 2 (n - 1) + 2p(Sn 2 - In + 2) - (2n - l)(4n - 1) . (3.7) 
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For n = 1 we get only the already excluded case p = 1/2. The n = 1 case 
(i.e. L = R 2 ) also has a de Sitter space-time solution, but since de Sitter is 
not self-similar, and since here we are focused on self-similar solution there 
are no solutions which meet our conditions for n — 1. 

Next, for n — 1/2 (i.e. L = R) we get Einstein's theory, which includes 
only the already excluded cases p = and p = 1/2. This too, is clear, as 
Einstein's theory is only of second order in contrast to the other cases which 
are of fourth order. 

The last special case is n = 1/4: Here we get unexpectedly only the 
already excluded cases p = and p — 1/2. 

Let us now give the general solution of eq. (3.7): The already discussed 
value p = 1/2 is a solution for every value n, though a spurious one for 
n < 1/2. 

Let us explain this in comparison with the field equation (2.1): In the 
notation used there, we have now L R = 2nR 2n ~ 1 . For n > 1/2 we have L R = 
and L = for every space-time fulfilling R = 0. Thus, the fourth-order field 
equation (2.1) possessing 10 independent components in the general case, now 
degenerates to one single second order scalar field equation, namely R = 0. 
Sometimes, that value of R where this happens, is called a critical value of 
the curvature scalar. For n < 1/2 we have, on the other hand, a positive 
power of R in the denominator, when the expression DLr is calculated, so 
space-times with R = should not be considered in that case. 

In addition to the p— 1/2 solution of eq. (3.7) there is a second solution 
which depends on n and which we will denote by p . This other solution p , 
which has already been found and discussed in [28] and [32], is given by 

(2n-l)(4ra-l) ,,31 /0 oX 

2 (n -i) --'-'"-r—i < 3 - 8 > 

One immediately sees that p = for n = 1/2 and n = 1/4, as discussed 
before. For n — 1, as expected, no value po exists. The only cases where 
Po = 1/2, i.e., that eq. (3.7) has the double solution 1/2, are n = and 
n = 5/8. 

For 1/4 < n < 1/2 and for n > 1 we have p < leading to a sudden 
singularity (Big rip), i.e., a(t) — > oo and simultaneously R(t) — > oo as t — > 
t < oo for a finite value t . 

In order to get an idea of the behaviour of p as a function of n, we look 
for its local maximum and minimum by evaluating the equation dpo/dn = 
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from eq. (3.8). This occurs at n = 1 =F = {0.3876; 1.6124}, where the 

(— ) is related to the local minimum and the (+) is the local maximum. These 
values of n correspond to related values of p as follows p = — 5 ± 2\/6 = 
{-0.101; -9.899}. 

Putting these results together shows that in the range 1/4 < n < 1/2 
we have —0.101 < p < and for n > 1 we have p < —9.899. Power-law 
inflation (PLI) takes place if p > 1. The equation p — 1 gives with eq. 
(3.8) the solutions n = (1 ± v^)/4 = {0.683; -0.183}. Therefore, PLI takes 
place for n < —0.183 and in the interval 0.683 < n < 1. Again this indicates 
(in the range po < 0) that we may encounter a Big rip singularity. 

Qualitatively one finds that the general solutions of eq. (3.6) have either 
self-similar solutions or de Sitter space-time solutions (which occurs for n = 
1) as asymptotic attractors. To illustrate this for the case n = 1 (i.e. for 
L = R 2 ) we recall that a closed form solution was already found in reference 
[24]. This solution (with the condition T^-l) reads 

ds 2 = - (l + T 3 ) ~ 2 dT 2 + (l + T 3 ) ~ 2/3 • (dx 2 + dy 2 + dz 2 ) . (3.9) 
In the limit T — > ±oo it is asymptotically equal to 

rfs 2 = _ T S dT 2 + T -2 . ^2 + dy 2 + ^ 

This is easily seen to be isometric to the Friedmann radiation model, i.e. to 
eq. (3.1) with a(t) = \/t for t — > 0. In the other limit, namely T — > — 1, 
it proves useful to write T = r — 1 with a small and positive or negative 
value of r. In leading order in r one then gets from eq. (3.9) the asymptotic 
behaviour 

ds 2 = -(3r)' 2 dr 2 + (3r)~ 2 / 3 • (dx 2 + dy 2 + dz 2 ) . 

With 3t = — In \t\ we can show this to be the exact de Sitter geometry, i.e. 
eq. (3.1) with a(t) = e*; here the region r — y corresponds to the region 
t — > oo. The exact solution eq. (3.9) confirms the general idea that the 
non-self-similar solutions of eq. (3.6) connect two different self-similar (or 
one self-similar and one de Sitter) solutions as asymptotic attractors. 6 

6 The de Sitter space-times represent a one-parameter set of space-times parametrized 
by the Hubble parameter. Thus the difference between the case n = 1, which has de 
Sitter space-time solutions, and the other values of n with n < 1 is that for n < 1 every 
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For n = 5/8, the exceptional case where eq. (3.7) has the double solution 
p — 1/2, (see also subsection 4.2. below for a discussion of this case), a 
non-self-similar exact solution was given in [28]. Its metric reads 

ds 2 = -a 2 ((In ay^da 2 - dx 2 - dy 2 - dz 2 ) . (3.10) 

Both for a — > and for a — » oo, the In ct-term in this metric can be replaced by 
1 to get the asymptotic behaviour. That means in both limits the asymptotic 
behaviour is 

ds 2 = —a 2 (da 2 — dx 2 — dy 2 — dz 2 ) 

which is again, as expected, the Friedmann radiation model, i.e. eq. (3.1) 
with a(t) = y/i. 

A further non-self-similar exact solution from [28] exists for the case n = 
1/4, (see also eq. (4.17) below for a further discussion of this case). Its metric 
is 

ds 2 = -a 2 (exp(a 6 )da 2 - dx 2 - dy 2 - dz 2 ) . (3.11) 
For the remaining values n, [28] gives the solution 

ds 2 = -a 2 ((a 6 + l) 2 ^ h da 2 - dx 2 - dy 2 - dz 2 ) (3.12) 

with / = -2 + 1/po and b = (8n-5)/(2n-l). For n -> 1/4 one has / ->■ oo, 
and for n — > 5/8 one has b — > 0, thus explaining the extra role of these two 
values of n. 

The metric in eq. (3.12) is not easily transformed to a synchronized time 
coordinate since finding the transformation between equation (3.12) and (3.1) 
is not usually possible to do in simple, closed form. 

We now count the degrees of freedom. For a general solution of a fourth- 
order field equation and one free unknown (in this case a(t)) one expects 
4 initial values, namely a(0) and the first three derivative of a(t) at t — 0. 
The remaining information is contained in the field equations. However the 
00-component of the field equation is a constraint, reducing the order by 
one. Thus the general solution is expected to have 3 free constants - one 
is a ^translation, the second is multiplication of t by some factor, the third 
is multiplication of a(t) by a constant factor. In this sense, for every fixed 

initially expanding solution approaches the PLI-exact solution as t — > oo. And for n = 1 
one can formulate the attractor property of the exact de Sitter inflation as follows: every 
initially expanding solution approaches one member of the one-parameter set of de Sitter 
space-times as t — > oo. 



11 



value of n, eq. (3.12) is a general solution even if no free parameter is 
explicitly shown. It is possible to compensate both the multiplicative time 
and the multiplicative space as follows: The multiplication of t by a non- 
vanishing factor can be compensated by a scale-transformation of the metric 
bringing solutions to solutions as the field equation is scale-invariant, and the 
multiplication of a by a constant factor can be compensated by multiplication 
of x, y and z by a constant factor. 7 

Now we test in which sense eq. (3.12) can be used in the case n — 1. 
For n — > 1 we get p — > ±oo, i.e. / = — 2 and 6 = 3. Then metric (3.12) 
transforms to 

ds 2 = _ a 2 ( (fl 3 + ^-4/3^2 _ dx 2 _ dy 2 _ ^ 

In principle, this metric should be isometric to metric (3.9). A closer inspec- 
tion shows the following: this metric with a > is isometric to (3.9) in the 
region T < — 1. 

As final observation we mention the fact, that the limit n — > ±oo is 
possible for metric (3.12). It leads via eq. (3.8), p = =Foo to / = — 2 and 
6 = 4, i.e. to 

ds 2 = -a 2 ((a 4 + l^da 2 - dx 2 - dy 2 - dz 2 ) . 

In synchronized time this equals 

ds 2 = -dt 2 + sinh(2t)(rfx 2 + dy 2 + dz 2 ) 

which is asymptotically de Sitter space-time as t — > oo. This metric repre- 
sents an exact solution of the equation 7' + 47 + 27 2 = which itself is the 
outcome of the limit n — > 00 applied to eq. (3.6) after it has been divided 
by -An 2 . 

3.2 Cosmological solutions in G n and G In G-gravity 

For L = G n , i.e., with a = and (3 = 1, eq. (3.5) reduces to 

= (n - 1) ((1 + 7 ) 2 - n(Y + 4 7 + 4 7 2 )) . (3.13) 

Obviously, this equation reduces to = if we put n — 1 in accordance with 
the fact, that G represents a divergence. 

7 This latter possibility is a consequence of the fact that the Euclidean 3-space is self- 
similar. Neither the closed nor the open Friedmann models share this property. 
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However, we may also do the following: replace (3* — ■ (n — 1) in eq. 

(3.5) for n ^ 1. Then we apply the limit n — > 1 to eq. (3.5) and get after 
multiplication by —(1 +7) the following non-trivial result: 

= 3a(l + 7) (67 + 3 7 2 + 27) - 2/3* (l - 2 7 - 3 7 2 - 7) . (3.14) 

If we now put a — and j3* — 1 and divide by —2, we get the equation 

= l-2 7 -37 2 -7', (3.15) 

which is just the relevant equation for the Lagrangian L = GlnG, see eq. 

(4.6) of [17]. Of course, if we insert L = GlnG directly into eq. (3.4), we 
also get eq. (3.15). 

We now make the analogous calculations as in the previous section. To 
begin we insert 7 = — 1/p into eq. (3.13) for n 7^ 1 and get 

= p 2 -2p(l-2n) + l-4n (3.16) 

which possesses two solutions p — 1 and p — 1 — An. These solutions coincide 
with analogous solutions in [5] - see eq. (30). There is a degenerate solution 
p = 1 in the irrelevant case n — 0. For n = 1 we can insert 7 = —1/p into 
eq. (3.15) and get solutions p — 1 and p = —3. One gets the same results if 
one puts n — 1 directly into the previously mentioned solutions of eq. (3.16). 
This property once more confirms the correctness of the claim that, in this 
context, the limit of G n as n — > 1 is equal to GlnG and not G. Also the 
singularity of GlnG at G = is a mild one. 

There is no value of n for which a de Sitter solution exists for L = G n . 
For every value n^O, there exist two different self-similar solutions: a(t) = t 
and a(t) = t 1_4n . All other solutions are not self-similar. For 1/4 there 
are three of them, in case of n — 1/4 only two of them. Every of these not 
self-similar solution possesses these two self-similar solutions as asymptotic 
limits. 

Now we go into further details. For n(n - 1) / we divide eq. (3.13) by 
— n{n — 1) and get 

= 7 ' + 47 + 47 2 - (l + 7 ) 2 /n. (3.17) 

In the limit n — > ±00 we get, similar as in the preceding subsection, a sensible 
equation, namely 

= 7 ' + 47 + 47 2 . (3.18) 
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Besides the solution 7 = — 1 valid for all values n, we have the additional 
constant solutions: (i) for eq. (3.18) we have 7 = which is a de Sitter 
solution; (ii) for eq. (3.17) with n 7^ 1/4 we have 7 = l/(4n — 1). This 
becomes directly obvious if we rewrite eq. (3.17) as 

= 7' + (1 +7) • {{An- 1)7- I) /n. (3.19) 

For n — 1/4, eq. (3.17) simplifies by the vanishing of its quadratic part in 7, 
and we get 

= 7' - 4 - 4 7 . (3.20) 

The solution of eq. (3.20) reads 7(7-) = c\ • e 4r — 1 with an arbitrary constant 
C\. The case c\ = is already covered above by a(t) = t. For c\ 7^ we 
can achieve C\ = ±1 by a suitable translation of r corresponding to a scale- 
transformation of the metric. Inserting a = e T we then get for 7 from eq. 
(3.2) 

7 = ±a 4 - 1 = h h~ 2 = a d/(d) 2 - 1. 

This leads to ±aa 3 — a/a which can be integrated once to c ± a 4 / A = ln(d) 
with an arbitrary constant c. Thus 

a = exp(c ± a 4 /4). 

Obviously, by a suitable rescaling of x, y and z we can get c = 0, so the only 
two non-self-similar solutions for the case n = 1/4 are 

d = exp(±a 4 /4). (3.21) 

The final integration is not possible in closed form, but the asymptotic be- 
haviour is obvious from eq. (3.21): For a — > we have d — > 1, so the PLI- 
solution a(t) = t is the typical behaviour near the singularity. For a — > 00 
and the upper sign, d becomes very large, so we have a Big rip singularity 
in the finite future. For the lower sign, however, d — > 0, the scale factor a{t) 
remains regular along the whole positive t-axis and approaches a(t) — > 00 as 
t — > 00. 

Let now n be different from and 1/4. Then p — 1 and p = 1 — An are 
both different from each other and different from zero. The related values of 
constant 7 = — 1/p are 7 = — 1 and 7 = l/{An — 1). 

First case: 8 n > 1/4, then l/{An — 1) is positive. The 3 non-self-similar 
solutions can be characterized by the magnitude of the values, which the 

8 For the special case n = 1 this can already be found in section 4 of [17]. 



14 



parameter 7 is covering: one contains 7 = 0, the second one contains very 
large positive values 7, and the third one contains values 7 < — 1. They 
behave as follows: That one containing 7 = has an initial singularity (a — > 
0) of a(t) = t-type 9 and approaches a(t) = |t| 1_4n for a — > 00 representing a 
Big rip singularity. That one with very large values 7 starts from an initial 
singularity (but not of a(t) = t-type) and also approaches a(t) = |t| 1_4n for 
a — > 00. The third one, containing the values 7 < — 1 has an initial singularity 
of a(t) = t-type and approaches another singularity for large values of a(t). 
In sum: the exact solution a(t) = a ■ |t| 1_4n represents a local attractor for 
a — > 00, and the limit n — > 1/4 is compatible with the result a{t) — > a 
deduced there, see eq. (3.21). 

The second case is < n < 1/4, which is equivalent to l/(4n — 1) < — 1. 
One solution starts from a(t) = t at t — > and has another future singularity, 
the second one also starts from a(t) = t but approaches a{t) = t 1 " 4 ™ for t — > 
00. The third one starts from another type of singularity but also approaches 
a(t) = t l ~ An for t — > 00. Result: The exact solution a(t) = ao-t 1 " 4 ™ represents 
a local attractor for t — > 00. 

The third case is n < which is equivalent to — 1 < l/(4n — 1) < 0. The 
intermediate calculations are different from the preceding, second case, but 
the end results for the behaviour of a(t) are surprisingly the same. 

3.3 Cosmological solutions in aR 2n + /3G n -gravity 

Here we look for solutions of eq. (3.5) with a/3 7^ and ti^l. Due to the 
fact, that a and (3 are free constants, there is a multitude of solutions, even 
if it is clear that essentially only the quotient a/fl plays a role. 

To make any head way to getting some idea of what are the interesting 
solutions for this case one should find some motivation for choosing specific 
values of a, (3, n. One possible idea is the following: one could use that values 
of n, where the power t 1_4n for the G n model coincides with the power p from 
eq. (3.8) from the in i? 2n -model. However, this takes place for n = 1/4 which 
implies p = which is one of the excluded cases. For the present we leave 
further detailed consideration of this case to possible future work if /when an 
argument can be found to selecting a restricted set of the parameters a, (3, n. 

9 This statement means that in a series expansion of a(t) 7 the linear term is the lead- 
ing one near t m 0, which represents the minimal requirement for a phase of power-law 
expansion. 



15 



3.4 Cosmological solutions in aR 2 + (3G In G-gravity 

Now we put L = aR 2 + f3G In G into eq. (3.4) and get 

= 3a(l + 7) (67 + 37 2 + 2 7 ') - 2/3 (l - 2 7 - 3 7 2 - 7') , (3.22) 

which is identical to eq. (3.14) with f3* replaced by (5. We look for solutions 
of eq. (3.22) with a/5 7^ 0. To this end, we insert 7 = — 1/p into eq. (3.22). 
Without loss of generality we assume a — 1/3 which transforms eq. (3.22) 
to 

= (1 - 1/p) (-6/p + 3/p 2 ) - 2(5 (l + 2/p - 3/p 2 ) . (3.23) 

The leading term in the limit p — > 00 of eq. (3.23) is = — 6/p — 2(5. If 
we insert (5 = —3/p into eq. (3.23), then we get qualitatively the following 
result: A very small negative value of the parameter (5 in the Lagrangian 
L = R 2 /3 + (5GlnG leads to the replacement of the exact de Sitter solution 
from L = R 2 , which is a local attractor, to a PLI exact solution a(t) = t p = 
t~ 3 /P which also represents a local attractor. This shows how one can modify 
the dynamics from de Sitter to power-law inflation by the addition of the 
G In G-term. 

We also can do this without any approximation as follows: We ask for 
what values of the parameter (5 in the Lagrangian L = R 2 /3 + (5G\nG 
the scale factor a(t) = t p represents an exact solution: This occurs when 
/3 = -3(2p-l)/(2p(p + 3)). 

4 Solutions in terms of elliptic functions 

In this section we show that for a certain subset of the general isotropic 
cosmologies in fourth-order gravity given above it is possible to integrate the 
gravitational field equation (3.5) to get a closed form, exact solution for the 
Friedmann scale factor a(t). It is found that for a broad range of these cases 
where we are able to fully integrate eq. (3.5) that the scale factor, a(t), 
is given by elliptic functions - hypergeometric functions 2 -^i (a, b; c; z) and 
incomplete Gamma function, F(a,z). We try, as far as possible, to give a 
physical interpretation of these solutions paying attention to the early and 
late time behavior. 
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4.1 Solutions for the 7 = const, case 

Before moving to the solution to the general eq. (3.5) we first examine the 
case when d^jdr = in eq. (3.3). This implies 7 = const, so that eq. (3.2) 
becomes 

where k is an integration constant. Next using this solution for h(t) in 
a(t)/a(t) = h(t) and integrating gives 

ln(a(t)) = In (k - t 7 )~ 1/7 ->• a(t) — (k — t 7 )~ 1/7 , (4.2) 

thus we see in this case a(t) is a power law of t i.e. a(t) oc t n . So that a(t) 
will increase with t and also to avoid having k — t'j become negative for some 
t one needs to require 7 < so that for large t we will have a(t) oc (ItI*) 1 ^' 7 '- 
Thus we find that for 7 = const. a(t) is a simple power-law of t. 

In order for 7 = const, to be a possible solution one needs to insure that 
there is a solution to eq. (3.6) such that A7 2 + B'-f + C = 0, where A, B, C are 
constants determined by n, a, (3. If there is no solution to A7 2 + B'-f + C = 
then the 7 = const, is not possible. If it happens that 7 = —2 or 7 = 
—3/2 then the expansion law will be that of radiation or pressure- less dust 
respectively. 

4.2 Exact solution for the case L = aR 2n 

In the case when the Lagrange density is given by L = aR 2n the field equation 
becomes eq. (3.6) which can be written as 

y = A 1 2 + B 1 + C (4.3) 

where the prime is differentiation with respect to r and the coefficients 
A,B,C are given in terms of n as 

4w-l 16n 2 - 14w + 4 2n - 2 

2n ; 2n(2n - 1) 5 n(2n - 1) { } 

The solution to eq. (4.3) is 

r d ^ 2 , /b + 2a 7 \ 

t = / — — — = — arctan 4.5 

J A^ + B^ + C D \ D ) y ' 



17 



where D = \f—B 2 + AAC . Inverting eq. (4.5) gives 



D 



/Dt\ B 
U j ~ 2A 



(4.6) 



Now taking into account the form of A, B, C from eq. (4.4) one can show 
that for all n that -B 2 + AAC < -0 so that D = y/—B 2 + AAC is imaginary 
(for n = 5/8 the equality holds and D — 0). Thus in eq. (4.6) we replace D 
with iDi with D 1 = \/ B 2 — AAC and taking into account that tan(iDi) = 
itanh(Di) we find that eq. (4.6) becomes 



7(t) = 7 tanh 



(4.7) 



In the case when n = 5/8 we have 7(7") = —B/2A = const, which was already 
discussed. In detail for n = 5/8 we find 7(7-) = —2 which implies a(t) oc t 1 / 2 
which is the expansion law for a relativistic gas with an equation of state 
P/P = 1/3. 

Next, we use this 7(7-) from eq. (4.7) to solve the following equation for 
h = a/a (where the overdot is differentiation with respect to t) 



7(7-) 



d 



dt \h 



dr d /l 



dt dr \h 



ldh 

h dr ' 



(4.8) 



where we have used h = dr/dt in arriving at the final result. Now we integrate 
eq. (4.8) for 7(7") from eq. (4.6) which yields 



Hh(r)) = - J \n 



cosh 



B 
2A 



T 



(4.9) 



or solving for h{r) gives 



h{r) 



cosh 



exp I ~- A r 



(4.10) 



Now using h(t) 
(4.10) becomes 



a 
a 



exp 



a(t)/a(t) for the l.h.s and r = ln[a(t)] in the r.h.s., eq. 
'£>i / Aln(a) 



In -iM 



2 y 

(Di+B)/2 



+ exp 



+ a ( " Dl+B)/2 ) _lM 



exp 



73 

2A 



ln(a) 



(4.11) 
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Finally we can integrate eq. (4.11) giving 

fdt = t= (^ l/A J da\a~ l (a^ +B ^ 2 + a^ +B ^ 2 ) 



l/A 



(4.12) 



2 {A-i)/a a , B _ D i 2AD 1 + B-D 1 n / \ * , 



B — Di w I 2ADr ' A ' 2AD X 

where k is an integration constant and 2 Fi(a,b;c; z) is the hypergeometric 
function. The form of the hypergeometric functions 2-^1(0, b; c; z) depends a 
lot on the values of the coefficient a, b, c and since from the last equality in 
eq. (4.11) these coefficients in our case can have a wide range depending on 
n the detailed analysis of this case may take some time. 

Thus we leave the detailed discussion of this case for a future work. Suffice 
it to say that the exact, closed form solution for the scale factor a(t) given 
above has a wide range of possible behaviors. 

We briefly discuss the case when A = to show that in this case the wide 
range of possible behaviors of a(t) implied by eq. (4.12) is greatly reduced 
to the point where one can make concrete statements about the resulting 
solution. The condition A = requires n — 1/4, B — 6, and C = 12. Using 
A = we find that upon integration eq. (4.3) yields 

f rf 7 ln(C + B 7 ) ln(12 + 67) ( . 

T = JWTC = B = 6^ (413) 

Inverting eq. (4.13) gives 

7 (r) = V - 2 (4.14) 
Next we integrate eq. (4.8) for 7(7-) from eq. (4.14) 



ln(/ i (r)) = ^exp(6r)-2r, (4.15) 



or solving for h{r) gives 



h(r) = exp ^— j exp(-2r) , (4.16) 

Now using h(t) = a(t)/a(t) for the l.h.s and r = ln[a(t)] in the r.h.s. of eq. 
(4.16) gives 

i-M^y^^Mx) (417) 
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Integrating eq. (4.17) gives, compare also eq. (3.11), 

where k is an integration constant and F[a, z] is the incomplete Gamma 
function. 

In figure 1 below we used eq. (4.18) to plot a vs. t. (i.e. we have 
graphically inverted this equation by plotting a as the dependent variable). 
From this we see the a(t) diverges at finite t (i.e. t = 1) which implies a Big 
rip. In order to avoid such a late time singularity we might exclude this case. 



Figure 1: Plot of a(t) over t with the integration constant taken as k = 1.5 
so that a > 0. One sees that a(t) diverges at finite t. This implies a late-time 
Big rip singularity for this cosmology. 



4.3 Calculations for G r 



For the case when L = G n the field equations takes again the form eq. (4.3) 
but with different coefficients, namely 

1 - An 



n 



2(1 - 2n) . 1 

t>G = , Og = — 

n n 



(4.19) 



For these values we have Dq 
7(t) in eq. (4.6) becomes 

7g(t) 



-B g + AA G C G = Ai thus the solution for 



2n , , rt , 1 - 2n 

- tanh (2r) 

1 v ' 1 - 



(4.20) 



An - I 1 - An 

Following the same steps as led to eq. (4.10) we now find 

h(r) = [cosh (2r)] n/(4n - 1} exp f- 2(1 ~ 2ra M . (4.21) 

v i - 4n y 

This in turn leads to an equation for a/ a similar to eq. (4.12) which can be 
integrated as before to give 

'l\t>/(i-*») 



J dt = t = 



ra/(l-4n) 



/da a" 1 (a'/» + ad"*)/")' 

4- 4 ,rb i i ; - a4 ( i ') + fc (422) 



1 \ n/(l-4n) , ]_ 

a(t) 2F1 
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where as before k is an integration constant and 2-^1(0, b; c; z) is the hyper- 
geometric function. 



Figure 2: Plot of a(t) over t for L = G n for n = 2. All values n > 1/4 
give the same general shape. One sees that a(t) expands quickly for large t. 



4.4 Calculations for L = G\n(G) 

In this case one again deals with eq. (4.3) but now with very simple coeffi- 
cients found by putting n — 1 in eq. (4.19). This yields 7' = — 37 2 — 27 + 1 
and implies 

A-cinG = —3 ; B GhlG = —2 ; C GhlG = 1 (4.23) 

and D G \ nG = i4. Thus we can start from eq. (4.21) and with the appropriate 
replacements A — > A G i nG , B — > B G \ nG , C — » C G \ nG , D\ — > iA arrive at the 
equivalent equation as eq. (4.11) 

Integrating eq. (4.22) gives 

j dt = t = 2 1 / 3 J da [a- 1 (a + a- 3 )~ V3 ] = 2 1 / 3 a(t) 2 F, Q, ^; | -a 4 (t)) + fc 

(4.25) 

where k is an integration constant and 2-^1(0, &; c; z) is the hypergeometric 
function. The function a{t) vs. t looks similar to the case when L = G n as 
given in figure 2. 



4.5 Calculations for L = aR 2n + (3G n 

In this case eq. (3.6) is not of the simply form of eq. (4.3) and thus the 
analysis will be different and much harder. For general n it will probably not 
be possible to get closed forms like in the previous cases. 
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5 Summary and conclusion 



In the introduction we had introduced the notions of scale-invariant and al- 
most scale-invariant Lagrangians. Let us here repeat the definitions in words 
instead of formulas as we did in section 1: A Lagrangian is scale- invariant 
if by a change of the length scale, the Lagrangian will be multiplied by an 
appropriate factor. A Lagrangian is almost scale-invariant if by a change of 
the length scale, the Lagrangian will be multiplied by an appropriate factor 
together with the addition of a total divergence not changing the field equa- 
tion. In our circumstances it turned out, that every almost scale-invariant 
Lagrangian can be written as limit of a sum of scale-invariant Lagrangians 
and total divergences. We conjecture, that this coincidence, which was for- 
merly already found in 2-dimensional gravity - see [84] for a more general 
discussion of this point, has a more general range of applicability. 

A related notion is the notion of a scale-invariant action: The action 
I = J L^—gd A x is called scale- invariant, if I does not change at all by a 
change of the scale. Further, / is called almost scale- invariant, if / only 
changes through the addition of a topological term by a change of the scale. 
For the Lagrangians discussed above, only L = aR 2 leads to a scale-invariant 
action, and only L = aR 2 + r yG hi G leads to an almost scale-invariant action. 
This is a strong argument for further detailed study of the gravitational 
Lagrangian 

L g = A + R + aR 2 + f3C ijkl C ijkl + 7 G In G (5.1) 

is being worth considered in more details than done up to now. Of course, 
the term /3Cij k iC ljkl - see [82] - by itself has a scale-invariant action, but we 
did not consider it in this paper, as it has no influence to the field equation 
within the Friedmann models. The terms A and R are added here since in 
the weak field limit such terms appear effectively. Of course the inclusion 
of matter will change the behaviour of the cosmological solutions discussed 
in this paper, but in the early stages of the Universe, with matter in the 
form of dust or radiation, the behaviour of the solutions above will only be 
marginally modified by the presence of the matter. The behaviour described 
in this paper will thus essentially correctly describe the dynamics. 

Finally, let us reformulate one of the key results of this work given in 
subsection 3.4 which is closely related to analogous calculations done for 
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higher dimensions in [87]: For p > 0, the Lagrangian 

3 p + 3 2p 

has the spatially flat Friedmann with scale factor a(t) = t p as exact vacuum 
solution. For large values p this is a local attractor solution and it represents 
a model for power-law inflation. 
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